Kandula, JASA-EL 


Spectral attenuation of sound in dilute suspensions 
with nonlinear particle relaxation 

Max Kandula 

ASRC Aerospace, NASA Kennedy Space Center, FL 32899, USA 
max. kandula- 1 @ksc. nasa.gov 

Abstract: Previous studies on the sound attenuation in particle-laden flows under Stokesian drag 
and conduction-controlled heat transfer have been extended to accommodate the nonlinear drag 
and heat transfer. It has been shown that for large particle-to-fluid density ratio, the particle 
Reynolds number bears a cubic relationship with cotj (where <yis the circular frequency and 

r^the Stokesian particle relaxation time). This dependence leads to the existence of a peak value 
in the linear absorption coefficient occurring at a finite value of cot j. Comparison of the 

predictions with the test data for the spectral attenuation of sound with water injection in a 
perfectly expanded supersonic air jet shows a satisfactory trend of the theory accounting for 
nonlinear particle relaxation processes. 

PACS numbers: 43.20.Hq, 43.50.Nm, 43.50.Gf 

1. Introduction 

Sound attenuation in fluids, representing the dissipation of acoustic energy from a sound wave, 
occurs through a number of physical processes involving molecular viscosity, thermal 
conductivity, and other dissipative or relaxation processes. M In all these absorption mechanisms, 
acoustic energy is converted into thermal energy, and other mechanisms deflect or scatter 
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acoustic energy. 1 Stokes 5 developed the first successful theory of sound absorption due to the 
effect of molecular viscosity of the fluid (internal friction). Helmholtz 6 theoretically investigated 
the effect of viscosity on the sound attenuation in a circular tube, considering the effect of friction 
near solid boundaries. Kirchhoff 7 theoretically treated sound attenuation due to the viscosity and 
heat conduction effects in an unbounded fluid medium and in narrow tubes. The so-called 
classical sound absorption in fluids includes Stokes’ viscous contribution and Kirchhoff s 
thermal conduction contribution to attenuation in unbounded fluids. 

When a fluid contains inhomogeneities such as suspended particles (solid particles, drops 
and bubbles) additional viscous and heat conduction losses occur in the immediate 
neighbourhood of the suspended particles. ‘‘ 3 ’ 8 A comprehensive review of the physics and 
scientific history of acoustic interactions with particulate mixtures is provide by Challis et al. 9 
Referring to Fig. 1, the acoustic intensity/ of a plane wave propagating through an absorbing 
medium is expressed by 

/ = v“" (1) 

where xis the distance traversed, Iq the intensity at x = 0 , and a, the intensity attenuation 
coefficient for the medium. The quantity a, depends on viscosity, thermal conductivity, and other 
factors such as molecular relaxation. 

Sound propagation in aerosols and fog has been studied experimentally and theoretically 
by several investigators. Sewell, 10 in his pioneering work, theoretically considered the attenuation 
of sound in a viscous medium containing suspended cylindrical or spherical particles (obstacles) 
with perfectly rough surfaces. With the aid of a scattering formulation (theory), Sewell predicted 
the attenuation of sound by rigid particles suspended in a gas, assuming that the particles are 
immovable. Epstein, 11 in a theoretical treatment of the attenuation of sound by spherical particles 
suspended in liquids or gases, derived and extended Sewell’s result by permitting the particles to 
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move. Epstein and Carhart 12 additionally considered heat conduction effects, and found that for 
fogs the effects of viscosity and heat conduction are both important, and approximately additive. 
Their formulation was extended by Allegra and Hawley 13 to include liquid-liquid as well as 
liquid-solid systems (with a different model for each case). 

The effect of transport processes on the attenuation and dispersion of sound in aerosols 
have also been reported by Soo, 14 accounting for nonstationary effects including the Basset 
history (effect of particle acceleration on viscous drag) and added mass terms. Basset history and 
added mass terms were included in an elegant coupled phase formulation by Harker and 
Temple. 15 Coupled phase effects were also treated by Evans and Attenborough 16,17 in an 
extension to the work of Harker and Temple 15 to incorporate thermal conduction. Nonstationary 
effects, which become important at low particle Reynolds numbers, are further investigated by 
Gumerov ,8 . Chow 19 considered the attenuation of sound in dilute emulsions and suspensions 
with viscous dissipation and thermal conduction, and additionally included the effect of surface 
tension, which was shown to be important in the presence of bubbles (and negligible for systems 
with solid particles or liquid droplets). 

Temkin and Dobbins, 20 in their classical work, theoretically considered particle 
attenuation and dispersion of sound in a manner which illustrates explicitly the relaxation 
character of the problem. For rigid particles, the linear droplet drag and heat transfer are 
respectively obtained from 


Fp - (i'p Wg ) 

(2a) 

Qp=lm! p k g [r p -T s ) 

(2b) 


which correspond to the zero droplet Reynolds number limit ( Re p -> 0 ), where 

Re^ — / 7 gjiVg —Up^dpl^ig (2c) 
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In the above, u is the velocity, T the temperature, p the density, p the dynamic viscosity, k the 
thermal conductivity, and d the particle diameter. The subscripts g and p stand respectively for 
the gas and the particle, and the primes denote fluctuations from the mean. The properties 
p g , p g , k g and p p refer to the mean values. According to Temkin and Dobbins, 20 Stokes linear 

drag law can be justified for 0 < cotj « 1 , provided that p g / p p « 1 and [cod p / 8Vg J « 1 . 

Recent data 21 on sound attenuation in supersonic air jets containing suspended water 
droplets reveal that the linear absorption coefficient displays a spectral peak. The particulate 
relaxation models for the sound attenuation are all based on Stokes drag (linear drag law) and 
pure conduction limit (linear heat transfer), and do not account or explain this attenuation 
behaviour. This article, primarily based on Ref. 22, attempts to investigate this attenuation 
behaviour by considering nonlinear drag and heat transfer laws applicable to relatively large- 
sized droplets. 

2. Analysis for sound attenuation 

The present analysis extends the work of Temkin and Dobbins 20 for dilute suspensions to 
accommodate the nonlinear drag and heat transfer laws, which become important at high particle 
Reynolds numbers Remand at high frequencies. Only sound attenuation is considered here, with 

sound dispersion excluded from consideration. 

2.1 Nonlinear particle drag and heat transfer 

The theory of Temkin and Dobbins 20 or dilute solutions is applicable at low mass 
concentrations C m =nqm p I pg « 1 , where «q refers to the mean number of particles per unit 

volume of mixture, and w^the mass of one particle. Without any loss of generality the 
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attenuation of sound for large particle Reynolds numbers with nonlinear particle relaxation may 
be expressed as follows: 


a = 


r c 0 a ^ 

yC m G); 


<° T d\ 


\ + co 2 r 2 d] 


Hr - 1 


'PP 


'Pg J 


COT t \ 


\ + Q) 2 T 2 \ 


(3a) 


In the above the quantity a refers the attenuation per unit frequency per unit mass fraction, a the 
linear absorption coefficient or the amplitude attenuation coefficient (imaginary part of the 
wavenumber), c p the specific heat, cgthe speed of sound in the gas phase, and ythe isentropic 

exponent (specific-heat ratio). Note that a = aJ2. 

The relaxation times and r t \ correspond to those under nonlinear drag conditions 
(generally representative of large-sized particles). They are related to the relaxation times r d and 


r, by the relations 



*d\ = *dV\^p\ 

r,i = W 2 ( R e p ) 

(3b) 

where 

(Re P )- Cd\ tC D , 

y/ 2 ( Re p > P r ) = Nu\ / Nu 

(3c) 


with Cq\ standing for the nonlinear drag coefficient, and Nu\fox the nonlinear heat transfer. In 
the above, 

T d= d lp P /( 18 /0 ( 4a > 


and 


T, = 


m p c pp ^* c ppdpPp | 

' 3 1 

(c \ 
L PP 

2nd p kg 12 fi g c pg 

. 2 ; 

< C PS J 


Pr r d 


(4b) 


where Pr = c pg ju g !k g stands for the Prandtl number of gas. Physically the quantity is a 
measure of the inability of the particles to follow (respond to) the fluctuations in the fluid 


motion. 3 Likewise, the quantity r„ is a measure of the thermal response of the particles to follow 


the fluctuations in the temperature of the fluid. 
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The drag coefficient and the Nusselt number in Eq. (3c) are defined by 


C D ~2F D l\p g ^y g 


Nu = h g d p /k g 


( 5 ) 


where F D is the particle drag force, and h g the gas-droplet convective heat transfer coefficient. A 
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good approximation (curve fit) for the drag coefficient C q\ is recommended as 


C D i = — + p= + 0.4 

Re P 


The first term on the RHS of Eq. (6) is the Stokesian drag coefficient defined by 


( 6 ) 


Cf) = 24 / Re p 


( 7 ) 


which is valid for particle Reynolds numbers less than about one. Eqs. (5) and (6) yield an 
expression for \f/\ as 






Re, 


+- 
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1 + 




: + 0.4 


( 8 ) 


With regard to the particle heat transfer, a good correlation for heat transfer (by 
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conduction and convection) is expressed by the well-known Ranz-Marshall correlation' 

Nu x =2 + 0.6ReJ 5 Pr 0 ' 33 (9) 

In the pure conduction limit, we have Aw = 2 , so that the function y/i 'm Eq. (3c) becomes 

y/ 2 =Nu { / Am = 1 + 0.3 Re J 5 Pr 0 ' 33 (10) 

It now remains to find a relation for the droplet Reynolds number Re pin terms of the density 
ratio and the particle relaxation time. 
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2.2 Droplet Reynolds number 

The determination of particle Reynolds number required in the evaluation of the functions y/\ and 
y / 2 in Eqs. (8) and (10) respectively is exceedingly complex. There exists relatively little 
information on the dependence of particle Reynolds number on the particle characteristics in two- 
phase flows. A study of turbulent diffusion of droplets in a gaseous medium 25,26 indicates a 
plausible relationship of the form 

Re p =f[p g lp p ,coT d ) (11) 

Since Eq. (3a) is applicable to large particle to fluid density ratios, we postulate here that the 
particle Reynolds number depends only on the particle relaxation time, and is independent of the 
particle to fluid density ratio: that is 

R e p =f{“> T d) ( 12 ) 

The physical basis for this supposition may be justified as follows. Particles with a smaller 
relaxation time are able to respond better to the fluctuations in the fluid motion, which leads to 
reduced particle relative velocity, and thus a reduced particle Reynolds number. 

In the present investigation, the following power law relation is proposed such that a peak 
in the linear absorption coefficient (as indicated by the measurements) is realized (an exponent of 
2.5 or greater for the relaxation time ensures this peak): 

Re p = c(cdt d ) 3 ( 13 ) 

The adjustable constant cis determined from a correlation of the theory with the test data. A 
value of c = 10 is found to be satisfactory based on the data of Norum 21 for water droplets in a 
supersonic air jet. 
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3. Results and comparison 

3. 1 Effect of Nonlinear Particle Relaxation 

The effect of nonlinear particle relaxation on the absorption coefficient per unit frequency is 
demonstrated in Fig. 2a for comparison with the theory of Temkin and Dobbins 20 for dilute 
concentrations. The results suggest that below cotj ~ 0.7 the nonlinear particle relaxation effects 

are seen to be unimportant. Beyond cor^ »0.7, the nonlinear theory departs from the linear 
theory. 

Fig. 2b shows a comparison of the predicted linear spectral attenuation of sound under 
nonlinear particle relaxation. The theoretical result by Temkin and Dobbins 20 for dilute 
concentrations is also presented for a comparison. The theory suggests that at high frequencies 
the linear absorption coefficient predicted by the nonlinear theory decreases with increasing 
frequency in accordance with experimental trend for large particle sizes. 

3.2 Comparison with experimental data 

A direct comparison of the present theory with the measured spectral attenuation by Norum 21 for 
water droplets in a perfectly expanded (jet exit pressure equals the ambient pressure; thus shock- 
free) supersonic jet is displayed in Fig. 3. The data correspond to hot supersonic jet of air from a 
convergent-divergent (CD) nozzle operation at a jet total temperature T t = 867 K, and a jet exit 

Mach number Mj= 1.45. The jet Mach number is defined as Mj = uj Icj, where the subscript 

j refers to the nozzle exit conditions. The mass flow rate (maximum considered) of water to that 
of the jet is about 0.85. The angle 0is measured from the jet inlet axis. The jet exit Reynolds 
number Re y) = Ujdj / Pgj is about 1.3xl0 6 , where d j is the jet exit diameter. At this condition, 

supersonic turbulent jet mixing noise dominates upstream noise radiation, and Mach wave 
radiation dominates the downstream noise radiation. In the data, water is injected at 45 deg. The 
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data include spectra measured at angles of 45 deg, 90 deg and 135 deg, thus highlighting 
directivity effects. 

From the measured ASPL (reduction in Sound Pressure Level) at a given frequency due to 
water injection, the linear absorption coefficient is deduced as 

« /cf max = ASPL/ASPL max (14) 

where the subscript refers to the peak spectral reduction. The comparisons suggest that the 
proposed theory based on the nonlinear particle relaxation processes satisfactorily describes the 
measurements for the spectral attenuation of sound, indicating a spectral peak. The inclusion of 
directionality effect (dependence on the angle of observation) on the spectral absorption is 
beyond the scope of the present work. 

It should be remarked, however, that the data on supersonic jet noise considered here for 
comparison are not directly pertinent to plane waves of sound for which the theory has been 
developed. Since noise from turbulent jets may be regarded as a superposition of plane waves of 
differing frequency, the absolute values of the absorption coefficient in the test data could be 
different from that expressed by the theory. In this connection, the reduction of turbulent jet 
mixing noise reduction by water injection was theoretically considered by the author in Ref. 27 in 
conjunction with the scaling laws proposed for jet noise 28 . Finally, the effect of droplet 
evaporation on the spectral attenuation of sound, as indicated by the theory of Marble et al. 29 is 
left out of account in the present analysis. 

4. Conclusion 

The theory proposed here for sound attenuation in dilute suspensions with nonlinear particle drag 
and heat transfer is shown to satisfactorily represent the test data for noise reduction with water 
droplets suspended in a supersonic jet. It is found that the nonlinear particle relaxation processes 
are primarily responsible for reduction in the linear absorption coefficient at high frequencies. 
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Figure Captions 

Fig. 1 Attenuation of a plane sound wave in a gas-droplet mixture. 

Fig. 2a (color online) Predicted absorption coefficient per unit frequency with nonlinear 
particle relaxation processes. 

Fig. 2b (color online) Predicted linear absorption coefficient with nonlinear particle relaxation 
processes. 

Fig. 3 (color online) Comparison of the predictions for the linear absorption coefficient 
with test data of Norum 21 . 
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Fig. 1 Attenuation of a plane sound wave in a gas-droplet mixture. 
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Fig. 2a Predicted absorption coefficient per unit frequency with nonlinear particle relaxation 
processes. 
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Fig. 2b Predicted linear absorption coefficient with nonlinear particle relaxation processes. 
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Fig. 3 Comparison of the predictions for the linear absorption coefficient with test data of 
Norum 21 . 
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